CCNY-HEP-99/4 
hep-ph/9905569 



Plasmon interactions in the quark-gluon plasma 



D. METAXAS and V.P. NAIR 



Physics Department 
City College of the City University of New York 
New York, NY 10031 
metaxas@ajanta. sci. ccny. cuny. edu, vpn@ajanta. sci. ccny. cuny. edu 



Abstract 



Yang-Mills theory at finite temperature is rewritten as a theory of plasmons which pro- 
vides a Hamiltonian framework for perturbation theory with resummation of hard thermal 
loops. 



1 Introduction 



It is by now a well-recognized fact that the so-called hard thermal loops (HTL's) play a 
very important role in Yang-Mills theories at high temperatures [Q. Hard thermal loops 
and the corresponding effective action have been studied rather extensively over the last 
few years [@, ||, |], ||. Hard thermal loops describe Debye screening, Landau damping, etc., 
and a resummation of perturbation theory including their effects is essential to avoiding a 
class of infrared divergences due to interactions of the electrostatic type. While this is clear 
conceptually, the actual implementation of a resummed perturbation theory has not been 
easy || [7|. Further, even in a resummed perturbation theory, there are still infrared problems 
which can be cured only by incorporating magnetic screening effects as well |5|. Such effects 
are expected to be of order g 2 T where g is the coupling constant and T is the temperature. 
In order to analyze magnetic screening effects, one must first incorporate hard thermal loop 
effects (which are of order gT ) and obtain an effective theory valid at lower scales, of order 
g 2 T. There have been many related approaches to this question in recent literature, based on 
kinetic equations, thermal Feynman diagrams, etc. [^, ID], [IT]. In this connection, it is worth 



emphasizing that magnetic screening involves the wave properties of the gauge fields. The 
magnetic mass may be considered as the (dynamically generated) mass gap of the effective 
three-dimensional theory of the zero Matsubara frequency modes; in terms of counting of 
powers of h this would seem to be a classical effect. However, the relevant classical theory is a 
wave theory, not a particle theory. In particular it is easy to see that, in the low energy limit 
where the zero Matsubara frequency approximation is valid, the single-particle wavefunctions 
have significant overlap, similar to what happens in Bose-Einstein condensation. Classical 
particle descriptions or corresponding kinetic equations will not suffice to generate a magnetic 
mass. One would thus like to have a way of incorporating the HTL-effects within a formalism 
which can systematically treat dynamical correlation effects as well as thermal corrections 
and preferably quantum corrections as well. The natural candidate for this would be an 
action formalism where one can add and subtract electric and magnetic mass terms and 
develop a resummed perturbation theory. While this is systematic, despite calculational 
complexity, the effective action for hard thermal loops involves nonlocality in time and some 
conceptual issues are better addressed in a Hamiltonian approach. Since the theory at 
finite temperature lacks manifest covariance anyway, there seems to be no serious drawback 
to a Hamiltonian analysis. In this paper we set up the basic framework of a resummed 
Hamiltonian analysis. 

The perturbative eigenmodes of the plasma of gluons at finite temperature are the plas- 
mons. In addition to the two transverse polarizations, there is also a longitudinal polariza- 
tion, which is physical mode at finite temperature, obeying the Gauss law. Our approach is 
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to rewrite the theory as a theory of such plasmons. We work out interactions of the plasmons 
to the quadratic order in the coupling constant. The Coulomb interaction between plasmons 
shows the screened behaviour, as expected. Perturbative calculations with this Hamiltonian 
will be a HTL-resummed perturbation theory. This sets up the basic framework. The next 
step is to use this to calculate corrections to HTL-effects, some of which are under way. 



2 Hamiltonian analysis 

As mentioned in the introduction, our approach will be to simplify the Hamiltonian analysis. 
The action for a non-abelian gauge theory, with the HTL-effective action added, is given by 

S=~Jf 2 + k J dQ[d 2 x T S wzw (G) + - J d /L xTr{G- 1 d^GA + 

-A-d+GG- 1 + A+G^A-G - A + A_)] (1) 

where Swzw{G) is the standard Wess-Zumino-Witten (WZW) action for G. G(x,Q) is a 
unitary matrix field depending on x^ and on the unit vector Q, Q ■ Q — Q 2 — 1, and obeying 
G\x,Q) = G(x, -Q). Also in (|), d± = \{d Q ± Q ■ 8) and A± = \{A Q ±Q ■ A). The 
G-dependent term of the action includes integration over the angles of Q. The parameter k 
is given in the lowest order calculation by (iV + ^Np)T 2 /6, where iV is the number of colors, 
Np is the number of quark flavors and T is the temperature. For our purpose, k can be 
considered free parameter. 

The Hamiltonian analysis of ([!]) has been given elsewhere |12[]. With the gauge choice of 
A = 0, the Hamiltonian is given by 



H 



f d 3 x E2+ o B2 + — ( d 3 xdn (4 + j 2 ) (2) 



where E% = d A^ and Bf = e^ k {djA a k + ^f abc AlA^) are the usual nonabelian electric and 
magnetic fields. The currents J± are given by 

J+ = ^d+gg- 1 = -it a J a + 

J_ = -—G- 1 D_G = -it a J a _ (3) 

47T 



and are related by J+(x, —Q) = J-(x, Q). 
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The equal-time commutator algebra which supplements the Hamiltonian (0) is given by 



{E?(x),A%y)} = -z<f«%<f(£-jfl 
[E«(x),J b + m = l ^-Q l 5 ab 5{x-y) 

[Jl(x,Q),J b + (y,Q')} = - l ^-(Q.VJ ab - f abc A c (x))5(x-y)5(QiQ') 

+if abc J c + {x 1 Q)5{x-y)5{Q,Q') (4) 

In addition to the Hamiltonian (g) and the operator algebra (|J), we must also impose the 
Gauss law as a constraint selecting the physical states |^), i.e., we must require that G a \ip) = 
0, where 

g a = (D ■ E) a + [dtt (J? + J a _) (5) 



Equations (§, |], |[) define the theory. The rest of this paper will be devoted to the analysis 
of these equations. 

The currents J± obey the generalized Kac-Moody algebra of (f|). The first step in our 
approach will be to introduce a canonical set of variables <j) a (x, Q), TL (x, Q) that depend on 
Q as well which obey the canonical algebra 

[r(^Q),U b (y,Q')} = ^ b 5(x-y)6(Q,Q') (6) 

with all other commutators equal to zero. 

We are interested in solving the theory to a certain order in the coupling constant or 
power of the structure constants f abc , and so we can solve (f|) as a series in f abc . The 
solution for J" is given by 

r + = rr - ^-q ■ v<p a - ^-q ■ A a 

87T 4tt 

K -f abc <f) b Q ■ V<f) c + -f abc <j) b U c 



48tt r ^ r T 
+ ^f abc f brs (Jl r (p s (p c + (p c 'tt r (p s ) + --- (7) 

It is easily verified that this solves the algebra @ to the quadratic order in f abc . This 
also agrees with the expression (|3|) for J" in a suitable parametrization of G(x, Q) in terms 
of (f) a (x,Q). J" is given by Q — > — Q in (0). Since the unitary matrix G(x, Q) obeys 
G^(x,Q) = G(x,—Q) we must have (p(x,—Q) = —<f>(x,Q) and H(x, — Q) = — U(x,Q). 
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(Strictly speaking, this condition has to be imposed as a weak condition; however, as the 
comments in the concluding section make clear, this will not affect the reduction of the 
Hamiltonian which follows.) By direct substitution from (|7]) and keeping in mind the above 
mentioned property, we find the Hamiltonian to quadratic order in f abc as 



n 



d 3 x 



2 k Jn 



An IQix Jn 



(Q ■ V# 



+ I / f abc (U + ■ A)y b Q ■ V0 C 



7T 



+ 7T I f abc f ars 
6k Jn 



+ 



K 



5767T Jn 
1 



~ ; jabc jars 

12 Jn 



IT C r IP) 

f abc f ars ((j) b Q- V0 c )(0 r Q- V0 S ) 

b U c (f) r Q ■ A 3 + <pQ ■ A 5 / IT) + 



(8) 



We must now consider the implementation of the Gauss law. This can also be done as a 
power series in f abc . Notice that the Gauss law operator Q a has the form 



ga 

Go 



Go + Gt + G a 2 + --- 

V • E a - — [ Q-V<p a 
47r Jn 



An 

f abc (A b ■ E c + 



(9) 



Qq is of zero order in f abc , is of first order, etc. We have also used (0) to simplify fl5|). 
Our strategy for the Gauss law will be to find a unitary transformation S = e F which has 
the property of transforming Q a to Qq. In other words we need 



s~'g a s w Ql 



(10) 



where the weak equality means "up to terms proportional to C/q." Such a strategy seems to 
have been known to many people before; it has been used recently for Yang-Mills theories at 



m 


13 


14 


t( 



formalism. We first define a quantity 

A a = V ■ A a + 2 



o (Q ■ p) 2 



Q ■ VIT 



(11) 
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A a may be considered as the gauge-fixing constraint conjugate to Qq. Indeed we find 

[gftx), A\y)] = + m 2 D )S{x - y) (12) 

where mr, = \/2k is the Debye mass. The inverse of the operator on the right hand side of 
(|12D will be denoted by G(x,y), i.e., 

f d 3 p e ip '( x ~ y ^ 

G(x,y) = T (13) 

J {2-k) 4 p 1 + mr D 

Defining 

X a (x)= fd 3 yG(y,x)A a (y) (14) 



we see that 

ig a (x),x b (y)] = -i5 ab 5(x-y) (15) 

Writing S = e F ~ 1 + F, to the first nontrivial order, (^) reduces to 

[g a (x),F]+g<?*o (16) 

The solution to this equation is given by 

F = -i J X a f abc (A b -E c + J^ b IT c ) - ~J x a diX h f ahc n + • ■ ■ 
Et = EZ-±f a QiP (17) 

(Notice that V • E = Qg.) 

The physical states can now be easily constructed. Let \ipo) be any state obeying the 
condition ^olV'o) = 0. Such states, as we shall see shortly, can be constructed by plasmon 
creation operators acting on the vacuum. Equation (|10D then shows that physical states \ip) 
may be written as 

\il>) = S|Vo> (18) 
In terms of such physical states we have for the matrix element of TC 

mn\^ 2 ) = {^ m \s- l ns\^ m ) (19) 

so that, in terms of the perturbatively constructed states \ipo), the Hamiltonian is effectively 
given by 

H eff = S~ l HS (20) 



6 



Qq and A a or (x a ) will commute with the creation-annihilation operators for the plasmons 
and, after obtaining TC e ff, we can set Qq = 0, A a = 0. Using the expression (Q) for the 
Hamiltonian and F as given by equation flT7|) we find 



with 



and 



where 



Ti 



eff 



Tin + Ti 



int. 



Tin 



Ti 



iiit 



E 2 



+-/(n-fg-A) 2 + -f /(g-v# 

K JO. 47T l07T Jn 

I f^diA^A^ + f abc f ars A b i AjAlA s j 
+~ / f abc (Tl + ■ A)y b Q • V0 C 

3 J 87T 

A f / a6c (n + ■ A) a \(f) b G(x, y)p c (y) 



; Q-v x G(x,y)p b (y)Q-v x cl> c (x 



(Q-p) 2 

+ 2L f f abc f ars (f) b U c (j) r U s 
3k J 

+ 57^ / f abc f ars ^ b Q ■ W c )(0 r Q ■ v< ^ s ) 

-—I f abc f ars {<p b Tl c <p r Q ■ A s + <j) r Q ■ A s <f) b Il c ) 
12 Jn 

' p a (x)G(x iy )p a (y) + --- 



f ab °(A h -E e + f ^ftjjcx 

Jn 



p a 



(21) 



(22) 



(23) 



(24) 



The interaction part of the Hamiltonian now contains the screened Coulomb interaction 
between plasmons. 
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3 Plasmon operators 

We now return to the construction of the plasmon operators which are eigenstates of the 
quadratic part of the Hamiltonian, viz., H . Such operators must evidently obey the equation 

[H ,a\ a (p)}=uj x (p)a{ a (p) (25) 

so that many-plasmon states are obtained by multiple applications of a\ a (p) on the vacuum 
obeying the condition a^(p)\0) = 0. In (|25"|), A is the polarization index. For a given spatial 
momentum p, we define a triad of unit vectors ef, A = 1, 2, 3, by 



c 3_ Pi 

Vp 1 

One can, if desired, make an explicit choice 

P2 'P. 



XX'a a 
6 e k 
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XX' 



(P3P1,P3P2,-(P1+PD) 



pI+pI JpI + pI 



{pI + pDp 2 

The plasmon creation operators may be taken to be of the form 



«i a (p) 



d 3 x e~ ip< 



e x ■ (uA a - iE° 



(An a + hr) 



(26) 



(27) 



(28) 



Substituting this in fl2"5j ) and simplifying, we see that solutions exist if uj\(p) are specific func- 
tions of p obeying certain dispersion relations. For the longitudinal plasmons, corresponding 
to ef , we find that ujl = uj^ia given by 

(Q • e 3 ) 2 



2^ 



dtt- 



For the transverse polarizations we find 



2 2 

uj t — p 



K 

2~k 



W£-(Q-p) 2 

(Q ■ e A ) 2 



<4 - (Q ■ pf 



(29) 



(30) 



The creation and annihilation operators, with the appropriate solutions for f\, f 2 substituted 



into fl28|) are, for A = 1, 2, 3, 



a x a (p) 
oclip) 



, 3 N x (p)e-^ x 



d a x 



V^xV 



, ■\x(l>)<'"-' 



d 6 x 



V^xV 



e x ■ (uA a - iE a ) + [dn / A (-uTT - ^(Q ■ p) 
x ■ (ooA a + iE a ) + [ dtt fx(-uoTl a + ^(Q- pfi 

J 87T 



2 xa\ 



(31) 
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where 



fx 



2 Q ■ e x 



Nx(p) 



u\-{Q- pf 



1-1/2 



(32) 



We have used plane- wave normalization appropriate to a cubical box of volume V = L 3 , so 
that p = (2ir/L)n, rii £ Z. The operators a,a^ obey the expected commutator algebra 



with all other commutators equal to zero. These operators create physical excitations in the 
sense that 



The operators (p)S are then operators which commute with Gauss law. We can also 

verify that A a (x) commutes with a, a 1 . Thus for matrix elements with states built up of 
a;t's, it is consistent to set Q§ = 0, A a = 0. 

We now turn to the question of expressing the Hamiltonian fl2"3|) in terms of the plasmon 
operators. The plasmons are collective mode excitations in the plasma. In general, one 
can write every operator, such as A, E, <p, II, in terms of a, plus operators /3, ft which 
represent other type of excitations. If the nonplasmon excitations are weakly coupled to 
a,a\ then an effective theory for the plasmons can be obtained by keeping just the a,a^- 
terms in A, E, (j),H. This is, in general, how we could obtain an effective theory for any 
kind of collective excitations. However, if a, form a complete set of operators, clearly 
this would not be an approximation, since there would be no independent /3, /3^-type modes 
anyway. As far as A^ Ei are concerned, it is easily seen that the plasmons do form a 
complete set. The operators 4>(x,Q), II (x, Q) have, in principle, an infinity of fields defined 
just on spacetime; these may be viewed as the coefficients of the expansion of <f>, U in terms 
of spherical harmonics constructed from Q. Thus completeness in terms of a, may seem 
doubtful. However, most of the degrees of freedom corresponding to the arbitrariness of the 
Q-dependence are irrelevant to us. This is easily seen by going back to the action ([J) and 
eliminating G by its equation of motion. As shown elsewhere, there are no new propagating 
degrees of freedom in G. In fact, the elimination of G leads to a purely A-dependent action 
and the analysis of this action shows that the only physical modes are the plasmons. (The 
action, with G eliminated, is nonlocal in time making a Hamiltonian analysis difficult; this 
is why we need the auxiliary field.) We may therefore use the plasmons as the only modes 




(33) 



[0o b (*W(p)] = o 



(34) 
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relevant to the interactions and reduce TC entirely in terms of a, cv . First of all, we write 

A1(x) as 

^W^^W+a^fp)] (35) 
Then, from (|3~3"D, we find a\{p) = — [a^ , Af(x)] and using fl31~|), 



We thus have 



axiip) 



ip-x 



In a similar way we find 

E?{x) = 

r(x) = 

n° = 



v Nx(p) 



'aS(p)^e-** + af(p)efe^ 
{-iu x ) \a a x (p)efe-^ - a?(p)e*e*- 



d 3 y K i {x,y)E^y) 



^ I d 3 y (Q-V) 2 K l (x,y)A«(y) 



where 



Ki(x,y,Q) 



d 3 p 



Jp-(x-y) 



E 



2Q ■ e A e A 



(36) 
(37) 

(38) 

(39) 
(40) 

(41) 



A Ul-(Q-P) 2 

Expressions ( |37j - ^0[ ) are easily checked to be compatible with the conditions Q$ = 0, A a = 0. 
The Hamiltonian (^) can now be expanded in terms of a, or in a completely straightforward 
way; we just have to substitute the mode expansions (]3"7|44"0"D in (123) . 



4 Discussion 

We shall conclude with some comments on the treatment of the auxiliary variables 0, II. 
The conditions Il(x, — Q) = — H(x, Q) and 0(x, — Q) = —(f>(x, Q) have to be imposed as 
weak conditions. Otherwise one may get inconsistencies in the application of the current 
commutation rules. Consider the expansion of the fields into harmonics on the two-sphere 
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given by Q. We may write 

— * _ — * 

(f}(x,Q) =2^4>i m (x)Y lm (Q) 

Im 

H{x,Q) ="£ll lrn {x)Y lm {Q) (42) 

l in 

The canonical commutation rules show that [0; m (x), Ui> m >(x')] = i5ii>8 mm >5(x — x'). The 
requirement of these fields being odd under Q — > —Q is equivalent to Ui m ps 0, <pi m ~ for 
even values of /. Ui m ~ for even / are first class constraints in the parlance of constraint 
analysis and the conditions <pi m ~ may be taken as the gauge-fixing constraints for them. 
One may set them to zero strongly one we have redefined commutators via Dirac brackets. 
Since the other fields x — (A, E i: (pi m , Hi m ) (for odd I) commute with the constraints, the 
Dirac brackets are the same as the Poisson brackets. Thus we may set the constraints to 
zero at this stage. 

In the Hamiltonian, the quadratic terms have no mixing of the even I- and odd l- 
components of the fields. In the interaction terms, there is mixing but we see that com- 
mutators of x's always generate terms involving the constraints and hence it is consistent to 
set the constraints to zero in the expression for the Hamiltonian. 

In summary, we have expressed the Hamiltonian with the HTL effects added in terms of 
plasmon creation and annihilation operators. One can use this for calculations of different 
types of higher order corrections; one can also use this as a starting point for variational 
calculations involving the thermal distribution of plasmons. Some of these are currently 
under way. 
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